Sensors 2008, 8, 8086-8103; DOI: 10.3390/s8 128086 



OPEN ACCESS 



sensors 

ISSN 1424-8220 

www.mdpi.com/journal/sensors 



Article 



Globally Optimal Multisensor Distributed Random Parameter 
Matrices Kalman Filtering Fusion with Applications 

Yingting Luo, Yunmin Zhu *, Dandan Luo, Jie Zhou, Enbin Song and Donghua Wang 

Department of Mathematics, Sichuan University, Chengdu, Sichuan, 610064, P. R. China; 
E-Mail: lyt83@163.com 

* Author to whom correspondence should be addressed; E-Mail: ymzhu@scu.edu.cn 

Received: 28 August 2008; in revised form: 26 November 2008 / Accepted: 3 December 2008 / 
Published: 8 December 2008 



Abstract: This paper proposes a new distributed Kalman filtering fusion with random state 
transition and measurement matrices, i.e., random parameter matrices Kalman filtering. It 
is proved that under a mild condition the fused state estimate is equivalent to the 
centralized Kalman filtering using all sensor measurements; therefore, it achieves the best 
performance. More importantly, this result can be applied to Kalman filtering with 
uncertain observations including the measurement with a false alarm probability as a 
special case, as well as, randomly variant dynamic systems with multiple models. 
Numerical examples are given which support our analysis and show significant 
performance loss of ignoring the randomness of the parameter matrices. 

Keywords: Random parameters matrices; Kalman filtering; Centralized fusion; Distributed 
fusion. 



1. Introduction 

Linear discrete time system with random state transition and observation matrices arise in many 
areas such as radar control, missile track estimation, satellite navigation, digital control of chemical 
processes, economic systems. Koning [1] gave the Linear Minimum Variance recursive estimation 
formulae for the linear discrete time dynamic system with random state transit and measurement 
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matrices without detailed rigorous derivation. Such system can be converted to a linear dynamic 
system with deterministic parameter matrices and state-dependent process and measurement noises. 
Therefore, the conditions of standard Kalman Filtering are violated and the recursive formulae in [1] 
can not be derived directly from the Kalman Filtering Theory. In this paper, a rigorous analysis 
(mainly in the appendix) shows that under mild conditions, the converted system still satisfies the 
conditions of standard Kalman Filtering; therefore, the recursive state estimation of this system is still 
of the form of a modified Kalman filtering. Reference [5] shows that this result can be applied to 
Kalman filtering with uncertain observations, as well as randomly variant dynamic systems with 
multiple models. 

Many advanced systems now make use of large number of sensors in practical applications ranging 
from aerospace and defense, robotics and automation systems, to the monitoring and control of a 
process generation plants. An important practical problem in the above systems is to find an optimal 
state estimator given the observations. 

When the processing center can receive all measurements from the local sensors in time, centralized 
Kalman filtering can be carried out, and the resulting state estimates are optimal in the Mean Square 
Error (MSE) sense. Unfortunately, due to limited communication bandwidth, or to increase 
survivability of the system in a poor environment, such as a war situation, every local sensor has to 
carry on Kalman filtering upon its own observations first for local requirement, and then transmits the 
processed data-local state estimate to a fusion center. Therefore, the fusion center now needs to fuse all 
received local estimates to yield a globally optimal state estimate. 

Under some regularity conditions, in particular, the assumption of independent sensor noises, an 
optimal Kalman filtering fusion was proposed in [11-12], which was proved to be equivalent to the 
centralized Kalman filtering using all sensor measurements; therefore, such fusion is globally optimal. 
Then, Song [7] proved that under a mild condition the fused state estimate is equivalent to the 
centralized Kalman filtering using all sensor measurements. 

In the multisensor random parameter matrices case, sometimes, even if the original sensor noises 
are mutually independent, the sensor noises of the converted system are still cross-correlated. Hence, 
such multisensor system seems not satisfying the conditions for the distributed Kalman filtering fusion 
given in [11-12]. In this paper, it was proved that when the sensor noises or the random measurement 
matrices of the original system are correlated across sensors, the sensor noises of the converted system 
are cross-correlated. Even if so, similarly with [7], centralized random parameter matrices Kalman 
filtering, where the fusion center can receive all sensor measurements, can still be expressed by a 
linear combination of the local estimates. Therefore, the performance of the distributed filtering fusion 
is the same as that of the centralized fusion under the assumption that the expectations of all sensor 
measurement matrices are of full row rank. Numerical examples are given which support our analysis 
and show significant performance loss of ignoring the randomness of the parameter matrices. 

The remainder of this paper is organized as follows. In Section 2, we present the concept of random 
parameter matrices Kalman filtering. In Section 3, we present an optimal Kalman filtering fusion with 
random parameter matrices and show that under a mild condition the fused state estimate is equivalent 
to the centralized Kalman filtering with all sensor measurements. In Section 4, we show that the result 
can be applied to Kalman filtering with uncertain observations as well as randomly variant dynamic 
systems with multiple models. More importantly, we will see that the Kalman filtering with false alarm 
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probability is a special case of Kalman with random parameter matrices. A simulation example is 
given in Section 5. And finally, in Section 6, we present our conclusions. 

2. Random Parameter Matrices Kalman Filtering 

Consider a discrete time dynamic system: 

x k+l = F k x k +v k (1) 
y k = H k x k +co k k = 0,1,2,... (2) 
where x k e R r is the system state, y k e R N is the measurement matrix, v k e R r is the process noise, 
and<y t e R N is the measurement noise. The subscript k is the time index. F k e R rxr and H k eR Nxr are 
random matrices. 

We assume the system has the following statistical properties: {F k ,H k ,v k ,co k ,k = 0,1,2,...) are all 

sequences of independent random variables temporally and across sequences as well as independent 
of*,,. Moreover, we assume x k and {F k ,H k ,k = 0,1,2,...} are mutually independent. The initial state x 0 , 

the noises v k , co k and the parameter matrices F k , H k have the following means and covariance. 

E(x 0 ) = ^ 0 ,E(x 0 -^ Q )(x () -^ 0 ) T =P 0 , (3) 

E(v k ) = 0,E(v k v T k ) = R Vt ,E(co k ) = 0,E(cD k c» T k ) = R ak (4) 

E(F k ) = F k ,Cov(f i ;,f* n ) = C fUL (5) 

£(^) = ^,C 0 v(^,^) = C AX (6) 

Where /* and ^ are the (i, j) th entries of matrices F k and//^ , respectively. 
Rewriting i^and H k as: 

F k =F k +F k (7) 

H k =H k +H k (8 ) 

And substituting (7), (8) into (1), (2) converts the original system to: 

x k+l = F k x k +v k (9) 

y k =H k x k +d) k (10) 

where 



G>k =( °k+H k x k 



(11) 



System (9), (10) has deterministic parameter matrices, but the process noise and observation noise 
are dependent on the state. Therefore, this would not satisfy the well-known assumptions of the 
standard Kalman filtering apparently. 

In the appendix, we give a detailed proof that system (9) and (10) satisfy all conditions of the 
standard Kalman filtering and derive the recursive state estimate of the new system as follows: 
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Theorem 1. The Linear Minimum Variance recursive state estimation of system (9), (10) is given by: 

= X k+\\k + Kk+l {y k+l~H K+l X k+l\k ) 



K k+1 - P k + Hk H Ll { H k + l P k+l\k H Ll +R a, k ) 



P k+i\k E k P klk F k + Ry t 



k+Uk 



P k+m+i K k+l H k+l jP k 
R, t =K t +E(F k E(x k x T k )F[) 
R 3t =R ak+ E(H k E(x k x T k )H T k ) 
E(x k+l x T k+l ) = F k E(x k x T k )F k +E(P k E(x k x T k )F[) + R Vt 

Xqk) — Ex 0 , P aQ — Var ( x 0 ) , E ^ x 0 x 0 j — Ex 0 Ex (j + P aQ 
where the superscript "+" denotes Moore-Penrose pseudo inverse, x k+m denotes the one-step 
prediction of x k+l , P k+m denotes the covariance of x k+m , x k+m+l denotes the update of x k+l and P k+M+l 



denotes the covariance of x k+m+1 . 



Compared with the standard Kalman filtering and noting the notations in (5), (6), the random 
parameter matrices Kalman filtering has one more recursion of E ( x k+l x T k+l ) as follows: 

E{x M x T k+l ) = F k E(x k x\ )F k +E(P k E(x k x T k )F k T ) + R Vt 

where 



and where is the (i, j) th entries of X k = E{x k x T k ) . 



3. Random Parameter Matrices Kalman Filtering with Multisensor Fusion 

In this section, a new distributed Kalman filtering fusion with random parameter matrices is 
proposed. The framework of the distributed tracking system is the same as those considered in [12-15]. 
The advantages of transmitting sensor estimates other than sensor measurements can be seen in [12- 
15]. We will show that under a mild condition the fused state estimate is equivalent to the centralized 
Kalman filtering using all sensor measurements. Therefore, it achieves the best performance. 

The / -sensor dynamic system is given by: 

Xk +l = F k x k+Vk> £=0,1,... 

y i k =H i k x k+ 4, i = l,..J (12) 
where x k e R r is the state, y[ e R N is the measurement matrix in i -th sensor, v k e R r is the process 
noise, and co\ e R N is the measurement noise. Parameter matrices F, and Hi are random. We assume 
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(13) 



that {F k ,H l k ,v k ,co J k ,k = 0,1,2,...} , i,j = \,...l, is a sequence of independent variables. Every single 

sensor satisfies the assumption in the last section. 

Convert system (12) to the following one with deterministic parameter matrices: 

x k + i = Fk x k+?k> A: = 0,1... 

y[=H[x k +a? k , i = l,...,J 
where 

fk= v k + Fk x k 

4 =G >i+Hi x k 

The stacked measurement equation is written as: 

y k =H k x k +d) k 

where 

and the covariance of the noise cb k is given by: 

Var{a k ) = R k 

Consider the covariance of the measurement noise of single sensor in new system. By the 
assumption above, we have: 

E (44 ) = E(4+H i k x k )(al+ H[x k f 

= E(44 +44H[ +Hix k co[ +H[x k x T k H[) 

= E(4co[) + E[H i k E(x k x T k )H['_ 

As shown in the last part of Section 2, every entry of the last matrix term of the above equation is a 
linear combination of E(hyh^) . Hence, when H' k and H[ are correlated, in general, 

E H l k E{x k x\ ^H ] k ^ 0 . Therefore, even if E [co l k co } k j = 0 , i.e., the original sensor noises are mutually 

independent, the sensor noises of the converted system are still cross-correlated, i.e., R k is non- 
diagonal block matrix. 

Luckily, when sensor noises are cross-correlated, in [7], it was proven that under a mild condition 
the fuse state estimate is equivalent to the centralized Kalman filtering using all sensor measurments. 

According to Theorem 1 and the Kalman filtering formulae given in [8-10], the local Kalman 
filtering at the i -th sensor is: 

x m = x \\k-\ +K 'k{y'k -Hk x m-\ ) = ( 7 - K kH'k) x k ] k -i + K ly'k 

jy i n; Tji r>i 

K k -K\k H k K k 

where 

K=Var(4) 

with covariance of filtering error given by: 
or 

P k \l=P k \l+HikXH[ (15) 
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where 



X k\k-\ ~ ^k X k-l\k-\ 



Pk\k - E( x k\k x k){ x m x k) 

p m-i = E( x k\k-\~ x k)( x k\k-\~ x k) 

We assume that the system has the following properties: the row dimensions of all sensor 
measurement matrices H[ to be less than or equal to the dimension of the state, and all H[ to be of 

full row rank. In many practical applications, this assumption is fulfilled very often. Thus, we 
knowH;(H;) + =/. 

According to [7] and Theorem 1, the centralized Kalman filtering with all sensor data is given by: 

x m ~ x k\k-\ + K k yy k ~Hk x k\k-i) 
= (l-K k H k )x k]k _ l+ K k y k 



(16) 



}T D 1 



K k =P klk H' k R k l (17) 
where, R k ] (*z) is the i -th column block of R k [ . The covariance of filtering error given by: 



p 1 =P 1 +H T k R k l H k 



(18) 



or 



= P t + H HL k " (*0^ ( H k ) H'BfHl 

i=\ 

P ak =(l-K k H k )P k]k _ x (19) 



where 

X k\k-\ = Fk X k-\)k-\ 

P k\k = E{ X k]k ~ X k ){ X k]k ~ X k) 

P k\k-\ = ^{ X k]k-1 ~ X k){ X k]k-\ ~ X k) 

Using (15) and (18), the estimation error covariance of the centralized Kalman filtering is given by 
using the estimation error covariance of all local filters: 

p := p ;l +H T k ±R- k \n)Ri(w;) + (p;; -£) m 

i=\ 

Using (17): 

K k y k =P k]k H T k Y 4 R k l {n)y[ 

' 1 , - v - (21) 

To express the centralized filtering x klk in terms of the local filtering, by (14) and (15), we have 

H f k R? y[ = pC Ky\ 



= p r 



^-(t-K^KJ (22) 



= P' x l -P' x l 

k]k k\k k\k-\ k]k-l 
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Thus, substituting (19), (21) and (22) into (16) yields: 

P l x = P l x +HlYR- k 1 (*i)Ri(H[) + (p r> x i -P'V ) (23) 

tit m Ht-i wc-i * t—i * V / K \ K J \ tit tit tit-i tit-i J V^~V 
<=l 

That is to say that the centralized filtering (23) and error matrix (20) are explicitly expressed in 
terms of the local filtering. Hence, the performance of the distributed random parameter matrices 
Kalman filtering fusion is the same as that of the centralized random parameter matrices fusion. 

Remark 1. In this paper it is assumed that all sensor observations are synchronous. In practice, this 
may be very rarely true. However, in the past 20 years, such assumption was used very often in the 
track fusion community (for example, see [7, 11-13] among others). One of reasons for this is that it 
is easy to extend the results for synchronous distributed Kalman filtering fusion to the corresponding 
asynchronous case by letting x = x' and P' = P' when the i th sensor does not receive its 

J J ° tit tit-i m tit-i 

observation y[ at time k to reduce the asynchronous distributed tracking system to be synchronous. 

Remark 2.The distributed systems here and in [7, 11-15] have a fusion center and the local sensor 
should transmit x' , x' , P' and P' to the fusion center. Our purpose in this paper is to derive a 

k\k m-i k\k tit-i 

globally optimal distributed fusion algorithm equivalent to the centralized Kalman filtering as the 
fusion center can receive all sensor observations. In this framework, the system has not only the global 
optimality, but also the good survivability in a poor situation. Clearly, such systems are different from 
the system considered [16]. The system in [16] does not require any form of central processing facility, 
and the algorithm there is highly resilient to loss of one or more sensing nodes, but costs more 
communication and has no real-time global optimality. Thus, both of them have their own advantages 
and disadvantages. The Random Kalman filtering in the framework of [16] may be another future 
research direction. 



4. Applications of Random Parameter Matrices Kalman Filtering 

In this section, we will see that the results in the last two sections can be applied to the Kalman 
filtering with uncertain observations as well as randomly variant dynamic systems with multiple 
models. 

4.1. Application to a General Uncertain Observation 

The Kalman filtering with uncertain observation attracted extensive attention [2-4]. There are two 
types of uncertain observations in practice. The first one is that the estimator can exactly know whether 
the observation fully or partially contains the signal to be estimated, or just contains noise alone (for 
example, see [2]). By directly using the optimal estimation theory, the Kalman filter for the first type 
of uncertain observations can be derived easily. The other uncertain observations belong to the second 
type, i.e., the estimator cannot know whether the observation fully or partially contains the signal to be 
estimated or just contains noise alone, but the occurrence probabilities of each case are known. 
Clearly, the latter is more practical. By applying the random measurement matrix Kalman filtering, we 
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can derive the Kalman filter with the second type of uncertain observations, which is much more 
general than that in [2-4]. 
Consider a system: 

x k + i = F k x k+ v k (24) 

^ZVo^+XV-'}^ (25) 

where all the parameter matrices are non-random and a set of multiple observation equations is 
selected to represent the possible observation case at each time. The random variable y k is either 

observable or unobservable. If y k = i, the measure matrix is H' k and the observation noise corresponds 

to o)' k . When the value of y k is observable at each time k , this is an uncertain observation of the first 

type and the state estimation with measurement equation (25) is converted to: 

y k =Hlx k +4 (26) 

which is obviously the classical Kalman filtering, i.e., the least mean square estimate using the various 
available observation of y k . To show the applications of the random measurement matrix Kalman 

filtering, we focus on the second type of uncertain observation, i.e., in (25), y k is unobservable at each 

timefc , but the probability of occurrence of every available measurement matrix is known. 

Consider that in (25), y k is unobservable at each timefc , but the probability of the occurrence of 

each measurement is known. Obviously, (2) is a more general form of (25) because only expectation 
and covariance of H k in (2) are known other than its distribution. The expectation of H k can be 

expressed as: 

H k =f.p J k Hi (27 ) 

H k =H' k -H k , with probability p[ (28) 

All that remains in order to apply the random measurement matrix Kalman filtering is just to 
calculate: 

= K k +E(H k E(x k x T k )H T k ) = R ak +f j p{{H[-H k )E{x k x T k ){H[-H k ) T (29 ) 

i=i 

Substituting (27}) and (29) into Theorem 1 can immediately obtain the random measurement 
matrix Kalman filtering of model (1) and (25). 

In the classical Kalman filtering problem, the observation is always assumed to contain the signal to 
be estimated. However, in practice, when the exterior interference is strong, i.e., total covariance of the 
measurement noise is large; the estimator will mistake the noise as the observation sometimes. In radar 
terminology, this is called a false alarm. Usually, the estimator cannot know whether this happens or 
not, only the probability of a false alarm is known. In the following, we will show that the Kalman 
Filtering with a false alarm probability is a special case of the uncertain observations of the above 
model (1), (25) are given. 

Consider a discrete dynamic process: 

x k + i = F k x k+ v k (30) 
y k =h k x k +a> k ,k =0,1,2,... ( 31 ) 
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where {F k ,h k ,v k ,co k ,k = 0,1,2,...} satisfy the assumptions of standard Kalman filtering. F t and h k are 
deterministic matrices. The false alarm probability of the observation is 1- p k . 

Then,we can rewrite the measurement equations as follows: 

y k =H k x k +(o k ,k = 0,1,2,... (32) 
where the observation matrix H k is a binary-valued random with: 

^{ H k=K} = Pk (33) 

Pr{H t =0} = l- ft (34) 

Due to (8): 

H k = PA (35) 

?r{H k ={l-p k )h k } = p k (36) 

Vx{H k =-p k h k } = \-p k (37) 

In the false alarm case, the state transition matrix is still deterministic, but the measurement matrix 
is random, by (35), (36) and (37), the covariance of the process and observation noises can be written 
as follows: 

R y t = R Vk (38) 
=K k +E(H k E(x k x T k )H T k ) = R 0}t +(l- p k ) p k h k E(x k x T k )h T k (39) 

Thus, the Kalman filtering with false alarm probability in this case is given by: 

X k+l]k+\ = X k+l]k + K k+1 ( y k+l — P k+ \h k +\ X k+\)k ) 
X k+l\k = F k X k)k 

K k+l = Pk +l p k+l )k h L (p 2 k + A +l p k +l \k h L + R a k ) + 

Pk+\]k = Fk^mFk +R v k 

P k+m+\ = (j ~ Pk+\K k+ \h k+ \)P k+ i\ k 

R„ t =K l +(i-Pk)pAE(x k x T k )h T k 

E ( X k + 1 X M ) = F k E ( X k X l) F k + K k 

x o]o ~ Ex o ' ^oio — V ar (-^o ) ' E (-^o^o ) — Ex 0 Ex 0 + P 0]Q 
In this section, we consider the application to a general uncertain observation for one sensor case. In 
a manner analogous to the derivation of Section 4.1, we can also give an application to a general 
uncertain observation for multisensor case using Section 3. The procedure is omitted here. 



4.2. Application to a Multi-Model Dynamic Process 



The multiple-model (MM) dynamic process has been considered by many researchers. Although the 
possible models considered in those papers are quite general and can depend on the state, but no 
optimal algorithm in the mean square error (MSE) sense was proposed in the past a few decades. On 
the other hand, when some of the MM systems satisfy the assumptions in this paper, they can be 
reduced to dynamic models with random transition matrix and thus the optimal real-time filter can be 
given directly according to the random transition matrix Kalman filtering proposed in Theorem 1. 
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Consider a system: 



x, 



■ k+l = F k x k +v k with probability p' k ,i = 1,...,/ (40) 

y k = H k x k+<°k (4i) 

where F k and v t are independent sequence, and H k is non-random. We use random matrix F k to 

stand for the state transition matrix. The expectation of F k can be expressed as: 

i 

F k=Hpl F k (42) 

7=1 

F k = F'l - F k , with probability p[ (43) 
A necessary step for implementing the random Kalman filtering is to calculate: 

R n = K +E(F k E(x k x T k )F k T ) = R Vk +f jP i k (F k i -F k )E(x k x T k )(F;-F k ) T (44) 

i=\ 

Thus, all the recursive formulas of random Kalman filtering can be given by: 

X k+l\k+l = X k+l\k + K k+i (v t+1 — FI k+] X k+l]k ) 
•"•t+lli = F k X k\k 

K M = P k+l\k H Ll ( H M P k + l\k H Ll +R co t ) 
P k+\\k = F k P k\k F k +R v k 

P k+\\k+\ = (j ~ K k+l FI k+l ) 

^ = ^ + Z Pi [ F l - ^ ) E [vl ) - F k J 
^oio — Fx {) , -^oio — V ar ( x o ) ' F (-^o-^o ) — Ex {) Ex Q + P m 



5. Numerical Example 



In this section, three simulations will be done for a dynamic system with random parameter 
matrices modeled as an object movement with process noise and measurement noise on the plane. The 
simulations give the special applications of results in the last section and show that fused random 
parameter matrices Kalman filtering algorithms can track the object satisfactorily. 

Remember that we have rigorously proved in Section 3 that the centralized algorithm using all 
sensor observations y[ at the fusion center can be equivalently converted to be distributed algorithm 

using all sensor estimates x l m . In addition, the computer simulations we have done show that the 

simulation results of two algorithms are exactly the same. It turns out that in the following numerical 
examples, we only compare the distributed random Kalman filtering and the corresponding standard 
Kalman filtering that ignores the randomness of the parameter matrices. Without loss of generality our 
examples assume the local sensors send updates each time when they receive a measurement. 

Firstly, we consider a three-sensor distributed Kalman filtering fusion problem with false alarm 
probabilities. 
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Example 1. The object dynamics and measurement equations are modeled as follows: 



x k + i= F k x k +v k 



y i k =H i k x k + 4,i = 1,2,3 

where [F k ,H\, v k ,co' k ,k = 0,1,2,...} satisfy the assumptions of standard Kalman filtering. The state 
transition matrix F, 



cos(2;r/300) sin(2^-/300) A 
-sin(2;r/300) cos(2;r/300) 



is a constant. The measurement matrix is given by: 



1 



,i = \,2,H\ = 



3^ 
3 



The false alarm probability of the / -th sensor is given by: 



1 



Pi 



0.01,1- pt =0.02,1- pi =0.03. 



The initial state x 0 = (50, 0) , P 0 ' l0 = / . The covariance of the noises are diagonal, given by 
R v = 1 , R , = 2 , i = 1,2,3. . Using a Monte-Carlo method of 50 runs, we can evaluate tracking 
performance of an algorithm by estimating the second moment of the tracking error, given by: 



1 

£2k= ^ 



I W _ II 2 



Figure 1 shows that the second moments of tracking error for three sensors Kalman filtering fusion 
without considering the false alarm (i.e. standard Kalman filtering) and three sensors random Kalman 
filtering fusion considering the false alarm (i.e. random Kalman filtering), respectively. It can be 
shown that even if the false alarm probability is very small, the distributed Random Kalman filtering 
fusion performs much better than the standard Kalman filtering. 



Figure 1. Comparison of standard Kalman filtering fusion and random Kalman filtering 
fusion. 

Est Err Var of Standard KF and Random KF 

12 



10 
8 



!i 



Random KF 

- - Standard KF 



T 



W 'i 



fi'di 1 2 
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In Example 1, both the sensor noises and the random measure matrices of the original system are 
mutually independent, so the sensor noise of the converted system are mutually independent. Now, we 
consider another example that both the noises and the random measure matrices of the original system 
are cross-correlated. 

Example 2. The object dynamics and measurement equations are modeled as follows: 



where i = 1, 2, 3 , { F k ,H l k ,v k ,a> k ,co' k , k = 0,1,2,...} satisfy the assumptions of standard Kalman filtering. 
The state transition matrix F k and the measurement matrices H[ are the same as Example 1 and co k is 

a large transition noise. When it happens, the sensors will mistake the transition noise as the 
observation. The false alarm probability of the transition noise is given by 1 - p k = 0.05 . Though the 

sensor noises co' k ,i = 1, 2, 3 are mutually independent and independent of m k , but the total measurement 

noises m\ are cross-correlated here. The covariance of the noises are diagonal, given 

by R v = 1 , R m = 2 R gi = 0.5 , i = 1, 2, 3. The initial state x 0 = (50, 0) , P 0 ' l0 = / . 

In this example, both the measurement noises and the random measure matrices of the original 
system are cross-correlated. Hence, the sensor noises of the converted system are cross-correlated. 
Figure 2 shows that the random Kalman filtering fusion given in Section 3 still works better than the 
standard Kalman filtering without considering the false alarm. This implies that the standard Kalman 
filtering incorrectly assumes that sensor noises are independent. 

Figure 2. Comparison of standard Kalman filtering fusion and random Kalman filtering 




y' k = H' k x k + co' k , with probability of p k 
= co k +a>' k , with probability of 1 - p k 



fusion. 
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Example 3. In this simulation, there are three dynamic models with the corresponding probabilities 
of occurrence available. The object dynamics and measurement matrix in (40) are given by: 
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F 2 = 



cos {in I 300) sin {in I 300) 

-sin (2^-/ 300) cos (2^/ 300) 

f cos (2W 250) sin {in 1 250) A 

-sin {In 1 250) cos (2^/ 250) 



with probability 0. 1 



with probability 0.2 



F 3 = 



cos [In 1 100) sin (2^-/ 50) 
-sin(2;r/50) cos(2;r/100) 



with probability 0.7 



1 1 



The covariance of the noises are diagonal, given by R v = 1 , R m = 1 . In the following, we compare 
our numerical results with the IMM. Since in this example, the occurrence probability of each model at 
every time k is known and mutually independent, it is also the transition probability in the IMM. 
Therefore, the transition probability matrix IT at each time in the IMM is fixed and given by: 

'0.1 0.2 0.7^ 
n= 0.1 0.2 0.7 
v 0.1 0.2 0.7 y 

Il(j, j) here means the transition probability of model i to model j . This assumption also implies that 

the model probability in the IMM is fixed as follows: 

,2-1=0.1, ^ 2 =0.2, ^ t 3 =0.7, k =1,2... 

Figure 3 shows that the random Kalman filtering given in section 4.2 still works better than the IMM 
with the fixed transition probability and model probability. This makes sense since the former is 
optimal in the MSE sense but the latter is not. However, in practice, the occurrence probability of 
each model is very often dependent on state or observation, and therefore not independent of each 
other in time. In this case, our new method offers no advantage. 
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6. Conclusions 

In the multisensor random parameter matrices case, it was proven in this paper that when the sensor 
noises, or the measurement matrices of the original system are correlated across sensors, the sensor 
noises of the converted system are cross-correlated. Hence, such multisensor system seems not to 
satisfy the conditions for the standard distributed Kalman filtering fusion. This paper propose a new 
distributed Kalman filtering fusion with random parameter matrices Kalman filtering and proves that 
under a mild condition the fused state estimate is equivalent to the centralized Kalman filtering using 
all sensor measurements, therefore, it achieves the best performance. More importantly, this result can 
be applied to Kalman filtering with uncertain observations as well as randomly variant dynamic 
systems with multiple models. The Kalman filtering with false alarm is a special case of Kalman 
filtering with uncertain observations. Numerical examples are given which support our analysis and 
show significant performance loss of ignoring the randomness of the parameter matrices. 



In the Appendix, we will give the proof of Theorem 1. Firstly, we give a detailed proof that system 
(9) and (10) satisfy all conditions of the standard Kalman filtering as follows: 

Lemma 1. Suppose random matrix F and random vector x are independent, then: 



Appendix 




Proof. By the properties of conditional expectation, we have that: 




Q.E.D. 



Lemma 2. 



{a)E{v k ) = 0,E{cb k ) = 0- 
(M)E(* 0 v[) = 0, (b2)E(x 0 a> T k ) = 0 
(cl)E(v k vf) = 0, (c2)E(a> k a>l) = 0, 
{d)E(v k v T k ) = R h ,E(a> k d> T k ) = R mt , 



(c3)E(v k d)f) = 0, Vk*l 



where 



\=K k +E(F k E(x k x T k )F k T ), 




Proof. 

(a): By the assumptions on the model (1) and (2) , and notations in (1 1) , it is obvious, 
(bl): Since {F k ,v k ,k = 0,1, 2,...} is independent of x 0 , 
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= E^x 0 x k F k ^ 

= E { x o{ F k-i x k-i+ v k-if F k} 
— E ^x Q x k _ l F k _ l F k ^ 

= E { X 0 { F k-2 X k-2 + ^-2 f F k-l F k } 

= E (^x 0 x k _ 2 F k _ 2 F k _ l F k ^ 

= Fix x T F T F T F T F T F T ) 
i^yx 0 x Q i Q i 1 ...i k _ 2 i k _ { i k j 

= E(xX) E ( F I) E ( F 1 T )... E ( F L) E ( F L) E ( F k) 

= 0. 

(b2): Similarly, E(x 0 a> T k ) = 0 . 

(cl): Without loss of generality, we consider the case of k > I only. 

E(v k vf) = E(v k vf + v k x]Pj + P k x k vf + P k x k x]Pj ) 

For x l is linearly dependent on F l _ l F l _ 2 ...F 1 F 0 x 0 ,v l _ l ,F l _ i F l _ 2 ...F l _ M v l _ i ,i = 2,3,...,/, and 

|F t v k ,k = 0,1,2,...} is independent of*,,, 

E(v k xfF l T ) = 0 
E(F k x k vf) = E{P k (F H x t _, +v t _ 1 )vf} 

= E ( F k F k-i x k-tf) 

= E { F k F k-l { F k-2 X k-2 +V k-l) V l} 

= E ( F k F k-l F k-2 X k-2 V f ) + E ( F k F k^k^f ) 

= E(F k F^..J?,xtf) + E(F k F^..J? M vtf^ 
= E(F k )E(F k _ 1 )...E(F M )E(v l vj) 
= 0. 

Noting that x l and {F k ,k = 0,1,2,...} are independent, by Lemma 1 we have: 

E(F k x k xjF l T ) = E{P k (F tAl +v t ,)^ ( r ) 

= E { F k F k-\ x k-i x i F i ) 

= E ( F k F k-l F k-2 X k-2 X f F ^ + F k F k-^k-2 xT l F l) 

= E(F k F k _ i F k _ 2 ...F l+l F l x 1 xfF l T ) + E(F k F^^^ 

= E { F k F k-l F k-2"- F l+l F l X l X l F l ') +E { F k F k-l F k-2— F l+l F l X l X l F l ) 

= E{E(f k F^ 1 ...F M )f l E( X rf)f?) 
= 0. 
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Hence, E (v k vf ) = 0 . 

(c2): Also consider the case of k > I only. Since x l is linearly dependent 
on F i-i F i-2~ F i F o x o> v i-i > F i-i F i-i- F i-M v ,-i > * = 2,3,...,/ , {F, , //, , v t , <y t , A: — 0, 1, 2, . . .} is independent of x 0 , 
and [F k ,H k , k = 0,1,2,...} is independent of ^.Moreover: 

) = E(co k cof + co k xjHf + H k x k cof + H k x k xfHf ) 

= E{H k {F k _ l x k _ l +v k _ 1 )xjHf} 
= F (H k F k _ 1 x k _ 1 xjHj) 

= E (Hk F k-i F k-2 x k-2 x fHf +H k F k _ 1 v k _ 2 xjHf} 

= F {H k F k ^ 2 ... F l+ ^x I x]Hj) + E(H 
= 0. 

(c3): 

E (v k ®i r ) = E [ V M + WJHJ + F k x k cof + F k x k xjHf ) . 

When k > I , 

E(v k a>f) = E(P k x k xjHf) 
= E{F k (F k _ l x k _ l + v k _ 1 )xfHf} 

= E ( F k F k-A-2 x k-2 x fHf ) + E(F k F k _ l v k _ 2 xfHf ) 

= E ( F k F k -A-2--- F M F i x i x fHf +F k F k _ 1 F k _ 2 ...F M v l xjHj) 
= 0. 

when k < I , 

E(v k xjHf ) = E[v k + v M ) r fff } 

= E {y k x L F i-iHf) 

= E (vk x l2 F * 2 Hl ) + E(v k vfi 2 F^Hf ) 

= E{v k x T k FlFl v ..Fl x H] ) + £ (vy k F^...F^Hf ) 
= 0, 
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E(F k x k xfHf) = E{P k x k (F w x w + tff } 
= E \h x k xT i-\ F i-iHj ) 

= £{F A (F A+ v,) r F/ +1 ...F^f} 

= E(F k x k x T k F k T F k T +v ..F [ T _ 1 Hf) 

= E(P k E(x k x T k )F k T E(F k T +1 ...F^Hf)) 
= 0 

) = + »v£ + W[ + F t ^F/ ) 

= K k +E(P k x k x T k F k T ) 

= R Vt +E(F k E(x k xl)F[) 

Similarly, f(^[) = ^ + E(H k E(x k x T k )H T k ). 

Q.E.D. 

By Lemma 2, system (9), (10) satisfies all conditions of the standard Kalman Filtering. Hence, we 
derive the recursive state estimate of the new system, i.e., Theorem 1 immediately. 
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